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A  summary  of^recent  results^  s  presented  on  the  subject  of 

,  A'' 

j, progressive  ply  cracking  in  fibrous  laminates.  -F+rst>>^valuation  of 
stiffness  changes  caused  by  systems  of  aligned  slit  cracks  which  are 
parallel  to  the  fiber  direction  in  a  unidirectional  composite  lamina  is 
discussed.  Results  obtained  by  the  self-consistent  method  are  presented. 
Next,  a  procedure  for  estimating  Instantaneous  crack  densities  and 
stiffness  changes  in  a  lamina  subjected  to  a  prescribed  strain  history 
is  outlined.  These  results  are  extended  to  analysis  of  laminated 
composite  plates  under  in-plane  stresses.  Specific  examples  and 
comparisons  of  analytical  and  experimental  results  are  presented  for 

two  graphite-epoxy  systems,  ±  (v  ‘A  W  $  it*  t 
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STATEMENT  OF  WORK 


In  many  fibrous  composite  systems  the  failure  strain  of  the  fiber 
far  exceeds  that  of  the  matrix.  Under  load,  the  difference  is  usually 
accommodated  by  matrix  cracking.  This  Is  frequently  observed  in 
monotonical ly  or  cyclically  loaded  laminated  plates,  where  each  layer 
may  contain  a  system  of  aligned  slit  cracks  which  grow  in  the  direction 
of  the  fibers  and  across  the  thickness  of  the  ply.  Such  cracks  are 
often  called  ply  cracks  or  transverse  cracks,  although  it  is  more 
appropriate  to  reserve  the  latter  for  cracks  which  are  perpendicular  to 
the  fiber  axis,  and  refer  to  cracks  which  grow  parallel  to  the  fiber 
direction  as  axial  cracks. 

In  polymer  matrix  composites  axial  matrix  cracking  typically 
starts  at  low  strain  levels  in  the  weakest  off-axis  ply.  As  loading 
continues,  cracks  appear  in  other  off-axis  plies;  also,  their  density 
increases  until  it  reaches  a  certain  saturation  level.  For  example,  in 
statically  loaded  (0/90)  graphite-epoxy  laminates  the  minimum  crack 
spacing  was  observed  to  be  equal  to  3. 5-4.0  ply  thicknesses  [1].  In 
metal  matrix  composites,  matrix  cracking  appears  to  be  caused  only  by 
cyclic  loads  which  exceed  the  shakedown  limit  of  the  laminate  [2]. 

Under  such  circumstances  the  matrix  experiences  cyclic  plastic  straining 
and,  consequently,  low-cycle  fatigue  failure.  Both  axial  and  transverse 
cracking  is  present,  the  former  in  off-axis  plies,  the  latter  in  zero 
degree  plies.  The  crack  patterns  and  densities  are  generally  similar  to 


those  found  in  polymer  matrix  systems.  However,  saturation  density 
increases  with  load  amplitude,  and  it  is  not  unusual  to  find  cracks  as 
close  as  one  ply  thickness. 

In  a  typical  part  of  a  laminated  composite  structure,  removed 
from  concentrated  loads  and  free  edges,  matrix  cracking  is  the  initial, 
low-stress  damage  mode  under  applied  load.  It  is  eventually  followed  by 
other  types  of  damage,  such  as  delamination  between  layers  and  fiber 
breaks;  but  these  appear  at  relatively  high  loads  which  may  exceed 
allowable  design  magnitudes.  In  contrast,  matrix  cracks  grow  at  low 
loads,  and  they  can  significantly  impair  stiffness  and  strength  of 
composite  structures,  especially  those  containing  many  off-axis  plies. 
For  example,  fatigue  tests  on  both  polymer  and  metal  matrix  laminated 
plates  indicate  that  stiffness  and  residual  static  strength  reductions 
caused  by  cracking  in  plies  may  be  as  high  as  10-50X  after  2  x  106 
cycles  of  loading  [2,3].  It  is  therefore  desirable  to  consider  the 
effect  of  matrix  cracking  on  composite  properties  in  design. 

Sufficiently  general  theoretical  models  of  progressive  cracking 
in  composite  laminates  are  apparently  not  to  be  found  in  the  literature. 
Such  results  as  are  available  for  angle-ply  laminates  have  been  obtained 
from  finite  element  calculations  [1],  while  other  studies  have  focused 
on  (Q/90n)s  laminates  [4], 

The  purpose  of  this  research  is, to  develop  a  procedure  for 


prediction  of  crack  densities  in  individual  plies  of  a  laminated 
composite  structure  as  a  function  of  applied  load,  and  to  evaluate  the 
effect  of  cracks  on  stiffness  of  the  structure. 

The  analysis  has  been  performed  in  the  following  steps: 

a)  Evaluation  of  overall  thermoelastic  properties  of  a  fibrous 
composite  which  contains  a  certain  density  of  aligned 

slit  cracks. 

b)  Evaluation  of  crack  densities  and  stiffness  changes 

in  a  single  ply  which  is  strained  In  a  prescribed  way. 

c)  Evaluation  of  crack  densities  and  stiffness  changes 
in  a  laminated  composite  plate  which  is  subjected 

to  prescribed  loading.  Specific  examples  were  solved 
for  laminated  plates  under  in-plane  loads. 

Comparisons  of  theoretical  predictions  was  made  with 
experimental  data. 


2. 


STIFFNESS  CHANGES  CAUSED  BY  A  SYSTEM  OF  CRACKS  IN  A  PLY 


2.1-  Preliminaries 


We  are  concerned  with  evaluation  of  overall  compliance,  thermal 
expansion  coefficients  (or  thermal  strain  and  stress  vectors),  and 
strain  and  stress  averages  in  the  phases  of  a  fibrous  lamina  which 
contains  aligned  slit  cracks,  and  is  subjected  to  uniform  mechanical 
loads  and  thermal  changes.  The  approach  to  the  problem  was  outlined  in 
reference  [5],  where  we  suggested  that  the  effect  of  matrix  crack 
systems  on  properties  of  fibrous  composites  can  be  analyzed,  in 
principle,  by  the  same  techniques  which  are  commonly  used  in  evaluation 
of  elastic  constants  of  composite  materials  and  fibrous  laminates,  e.g., 
by  the  self-consistent  method. 

The  essential  approximation  in  the  evaluation  of  stiffness  and 
compliance  changes  of  laminates  consists  in  the  replacement  of  a  cracked 
layer,  Figure  lb,  by  an  effective  medium  which  contains  many  cracks. 
Figure  la.  The  crack  densities  can  be  exactly  matched  to  provide 
identical  stiffnesses.  However,  the  cracks  in  the  layer  are  not 
entirely  surrounded  by  the  layer  material,  instead  they  interact  with 
neighboring  layers  which  have  different  elastic  properties.  This 
interaction  is  limited  to  the  vicinity  of  the  crack  tip,  thus  it  may  be 
important  in  analysis  of  local  crack  growth  at  the  interface,  but  it  has 


only  a  minor  effect  on  lamina  stiffness.  We  note  that  a  similar 
approximation  is  commonly  accepted  in  evaluation  of  elastic  moduli  of 
laminated  composite  structures  reinforced  by  monolayers  of  large 
diameter  fibers. 

When  the  composite  is  reinforced  by  monolayers  of  fibers,  such 
as  boron  or  silicon  carbide,  the  cracks  and  fibers  may  be  of  similar 
size.  The  appropriate  model  is  shown  in  Figure  2.  It  is  analogous  to 
that  of  Figure  1,  except  that  it  contains  three  phases  (fiber,  matrix, 
and  cracks),  whereas  the  model  of  Figure  1  can  be  reduced  to  two  phases 
(composite  "matrix",  and  cracks).  Accordingly,  the  models  in  Figures  1 
and  2  are  referred  to  as  two-phase  and  three-phase  models,  respectively. 

From  the  practical  standpoint,  the  effect  of  model  choice  on  composite 
stiffness  is  small.  The  simpler  two-phase  model  is  thus  sufficient  for 
analysis  of  all  fibrous  systems. 

In  the  three  phase  model  we  designated  the  fiber  as  phase  1, 
matrix  as  phase  2,  and  cracks  as  phase  3.  In  the  two  phase  model  phases 
1  and  2  are  joined  in  a  homogeneous  "matrix"  and  designated  as  phase  2; 
the  voids  or  cracks  remain  phase  3. 

The  self-consistent  analysis  of  the  cracked  composite,  as  outlied 
in  [5-7],  starts  with  a  composite  geometry  in  which  the  cracks  are  initially 
regarded  as  ellipsoidal  cylindrical  inclusions.  A  self-consistent 
estimate  of  overall  stiffness  of  this  medium  is  obtained.  Next,  the 


inclusions  are  evacuated,  i.e.,  replaced  by  voids.  Finally,  the  aspect 
ratio  of  the  voids  is  adjusted  so  that  in  the  limit  the  voids  change  to 
cracks. 

To  evaluate  crack  density,  the  cracks  in  Figure  1  are  first  replaced  by 
elliptic  cylindrical  voids,  with  a,  b,  denoting  the  major  and  minor  semiaxes. 
If  n  is  the  number  of  voids  per  unit  area  in  the  x-^-Dlane,  then  the 
volume  fraction  of  voids  is  equal  to  c^  »  iTabn,  and  c2  +•  c3  *  1.  Next, 
let  the  voids  be  reduced  to  cracks,  i.e.,  5  »  b/a  ♦  0.  Then 

C3  *  Tra2n<5  «  1  nQ6  (1) 

where  S  =  4na2  is  the  crack  density  parameter.  In  fact,  8  Is  equal  to 
the  number  of  cracks  of  fixed  length,  e.g.,  2a,  in  a  square  of  side  2a. 

For  example,  if  the  cracks  are  located  in  a  ply.  Figure  lb,  then  3 
measures  the  distance  between  regularly  spaced  cracks  in  terms  of  ply 
thickness  2a.  At  8  *  l,  the  distance  between  cracks  is  equal  to  2a,  as  3 
decreases  the  distance  increases  and  at  8  s  0  the  cracks  vanish.  We 
recall  that  the  observed  minimum  distances  between  cracks  in  a 
saturation  state  are  equal  to  3. 5-4.0  ply  thicknesses,  i.e.,  7a-8a  [1]. 
Corresponding  values  of  8  are  0.28-0.25,  but  values  as  high  as  3  *  1 
were  observed  in  the  8-A^  system  [2].  Therefore  0  <  S  <  1  is  the 
appropriate  range  of  3. 


On  the  macroscale,  the  cracked  unidirectional  composite  of 
Figure  1  can  be  regarded  as  an  orthotropic  homogeneous  solid.  The 
elastic  properties  of  the  "matrix"  are  identical  with  those  of  the 
fibrous  composite  and  can  be  easily  evaluated.  When  cracks  are 
introduced,  the  macroscopic  or  overall  elastic  moduli  of  the  solid 
change.  To  make  the  concept  of  overall  moduli  meaningful,  it  is 
necessary  to  consider  overall  uniform  loading.  Thus,  we  introduce 
uniform  overall  average  stresses  g  and  strainsT,  with  components 
arranged  in  (6x1)  column  vectors  and  related  by  constitutive  equations* 

a  *  L  s  ,  c  *  M  a  *  '2) 

where  L,  and  M  are  the  overall  stiffness,  and  compliance  (6x6)  matrices 
of  the  cracked  composite,  respectively.  M  *  L"1  when  the  inverse 
exists.  Effective  properties  of  uncracked  fibrous  material  (Ph*se  2)  are 
denoted  by  1°,  M°,  or  by  l^, 

*  As  in  [5-7],  (6x6)  matrices  are  denoted  by  capital  Latin  or  Greek 
letters,  e.g.,  L,  M,  A,  A,  8,  P,  Q,  and  (6x1)  matrices  by  lower  case 
bold  face  letters,  e.g.,  a,  e,  m,  £. 


For  a  dilute  concentration  of  aligned  slit  cracks,  the  matrices  L 
and  M  can  be  evaluated  by  the  self-consistent  method.  This  was  done  in 


Reference  [5];  the  results  are: 

L’L2’®"L2aL 
M  «  M2  ♦  7  a 

where  *  1°,  *  M°,  and 

8  *  J  *8 


(3) 

(4) 

(5) 


The  matrix  A  has  only  three  nonzero  components,  which  are 
expressed  in  terms  of  compliances  M  of  the  effective  medium  as: 


_  ^22^33*^23  ,  h  , 

*22  '  VLl  (al  V 


33 


A44  =  (M44  ^5^ 


(6) 


(M, 


66 


22M33’M23^(MnM33-M13)"  (fl£  +  ^ 


,2 


33 


where  and  are  the  roots  of 


r 


These  results  imply  that  only  three  compliance  coefficients 
and  Wg6  are  affected  by  the  introduction  of  cracks,  the  remaining  six 
terms  in  M  are  unchanged,  i.e.,  they  remain  equal  to  those  of  the 
uncracked  fiber  composite.  In  particular 


H11  *  •  m12'MU  ■  "13  "$3 


M23  *  **23  ’  M33  ’  ,M33  ’  M55  '  *55 


«22  -  ^2  +  8(M22M33  -  ^3><“?  +  “2>/M3 


H44  ’  ^  + 


M66  *  ^6  +  S C M22M33-Hl3 > ^ M11 M33 ^ 3 ^ 4  +  °2^M33 


and  from  (7): 


«1«2  *  (M11H33^3>/'M22M33-m223> 

“3  *  «2  *  CM33M66  *  2!M12M33-M13M33):l/CM22M33-H233 


The  unknown  shear  compliance  M44  can  be  obtained  directly  from 


••  .v. 


\  N  •«  \  \  \  ATVA' 


c 


M44  *  ^4  +  &  M55  *  (®2  ”55  *  4MW  M55)  ^ 


(U) 


The  remaining  unknowns  and  Mgg  are  found  from  (8)  and  (10). 

These  results  can  be  utilized  in  (4)  for  a  more  direct  evaluation 
of  the  three  nonzero  components  of  the  matrix  A  in  (6): 


r 

■ 


So  *  (Mfj  •  m  *  22  ,  44,  M)  (i2> 

Of  course,  the  same  result  follows  from  (6)  and  (7).  In  any 
case,  the  A  may  now  be  substituted  into  (3),  and  the  components  of  the 

J 

overall  stiffness  l  can  be  found  in  a  closed  form.  The  resulting 
expressions  are  given  in  reference  [7,  eqn.  (A-6)]. 


It  is  seen  that  the  results  are  remarkably  simple,  and  similar  to 
those  that  are  routinely  used  for  evaluation  of  elastic  moduli  of  an 
uncracked  composite  medium.  This  similarity  is  particularly  useful  in 
applications  of  the  results  to  laminated  structures. 


L. 
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Let  6  denote  a  uniform  temperature  change  from  a  reference 
temperature  90,  applied  to  the  composite.  The  constitutive  equation  (2) 
must  be  replaced  by 

<T  *  L  e  -  94  7  *  M  a  +  9n  ,  (13) 

'X.  %  %  % 

where  m  is  the  thermal  expansion  vector  and  4  the  thermal  stress  vector. 
From  (13)  one  obtains 


t  =  L  m  (14) 

It  is  probably  obvious  that  the  presence  of  cracks  does  not  affect 
free  thermal  expansion  of  the  composite.  Thus 


m 


(15) 


and  from  (3),  (4)  and  (14): 


l  =  (I  -  0UV)4O 

Explicit  forms  of  m  and  l  are: 


(16) 


-12- 


i 


(LU*  L12>  aT  *  l13  aA 


(L12+  L22}  °T  *  L23  °A 


<L13*  L23>  aT  +  l33  aA 


where  and  aA  are  linear  thermal  expansion  coefficients  of  the 
uncracked  composite  in  the  transverse  and  axial  directions,  respectively. 
Since  the  coefficients  of  1^  change  with  3,  £  Is  also  a  function  of  m, 
while  m  remains  constant. 


Again,  the  derived  expressions  are  valid  only  if  the  cracks  are 
open.  Closed  cracks  do  not  affect  thermal  response  of  the  composite. 


•iV'.  •  *  .  *  •  _ 


V*l  v*  O  O*  o  O  «>V%\ O O.  O  .  O.  KTr 


2.3  Local  Strain  and  Stress  Averages 


Certain  applications  of  the  above  results  require  information  about 
local  stress  and  strain  fields  in  the  phases.  Estimates  of  local  fields 
are  also  needed  to  distinguish  between  loading  conditions  which  lead  to 
either  open  or  closed  cracks  and  thus  delineate  limits  of  validity  of 
the  theory. 

It  is  clear  that  when  a  fiber  reinforced  composite  containing 
cavities  is  subjected  to  general  uniform  mecahnical  loading,  part  of  the 
applied  strain  (i.e.,  the  average  strain)  is  accommodated  by  cavity  ex¬ 
pansion.  Whereas  the  strain  of  a  typical  cavity  wall  becomes  unbounded 
in  the  limit  of  vanishing  aspect  ratio,  the  important  quantity  in  the 
macroscopic  study  of  composites  is  the  wall  strain  multiplied  by  the 
area  of  the  cross-section  (or  aspect  ratio)  which  tends  to  a  finite 
limit. 


From  the  formulae  given  in  [5,7],  and  in  Section  2.2  above,  it  is 
not  difficult  to  show  that  in  the  limit  of  slit  cracks  (5+0),  the 
overall  strain  7  is  given  by 


e  = 


eo 


+  c 


c  » 


(18) 
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where  eQ  is  the  average  strain  in  the  composite  “matrix"  : 

*o  *  *  *  A  L)g  +  8  A  L  m  8  .  (19) 

In  addition,  7C  is  the  crack  accommodation  strain 

tc  «  8  A  L(7  .  me)  ;  (20) 

it  represents  that  part  of  the  applied  (average)  strain  which  is 
accommodated  by  cracks. 

In  practice  equations  (19)  and  (20)  provide  us  with  estimates  of 
the  local  stress  and  strain  fields  within  the  composite. 

Viewed  from  another  standpoint,  equations  (19)  and  (20)  imnediately 
furnish  the  strain  concentration  factors.  As  far  as  the  matrix  is 
concerned, 

c  *  A  c  -  9a  (21) 

-o 

where 

A  *  I  -  sal  (22) 

-  *  (A0  ■  I)  *  •  8  A  l  m 


-15- 


It  is  perhaps  important  to  emphasize  the  different  physical 
interpretations  of  eQ  and7c*  Whereas  tQ  is  the  average  strain  in  the 
matrix,  t '  represents  that  part  of  the  overall  strain  which  is 

V* 

accommodated  by  all  the  cracks.  Therefore,  one  can  write  (20)  as 

?c'ic5'95c  •  (23> 


where 


A  »  8  A  L 
c 

a  «  FA  L  m 
~c  - 


(24) 


For  example,  it  is  immediately  obvious  that  the  average  thermal 
strains  in  a  fully  constrained  cracked  composite  ( e  *  0)  are 

lo  •  -  !c  •  7  A  L  m  6  . 

The  matrix  expands  at  the  expense  of  crack  closing. 


For  computational  purposes  It  is  advantageous  to  rewrite  ( 22-j ) 
with  the  help  of  (3)  in  the  form 

A  *  MQ  l  .  (25) 


It  is  often  convenient  to  have  explicit  forms  of  the  above  results 


The  nonzero  components  of  A  in  (22)  are: 


An 

**33  * 

A55 

*  1 

a21 

'  (M22  - 

m22) 

L12 

^22 

*  <«22  * 

M22) 

42 

^23 

*  <■$** 

m22) 

L23 

A44 

*  ">* 

A66 

Again,  are  compliances  of  the  uncracked  composite. 


(26) 


The  components  of  thermal  strain  concentration  factor  a  in  (22)  are 
all  zero,  except  for 


a22  "  "  ®  A22  ^L12  +  L22)oT  +  L23  °A^  <27> 

The  average  strain  ec  (20),  (23)  in  the  cracks  has  the  following 
nonvanishing  components: 


*22  *  *  A22CL12(eir*Te)  *  L22(c22'V)  +  L23(T33*aAe)^  (28) 

T23  *  (1  *  M44/M44^c23  r12  "  ~  M?6/,W*12 

We  note  that  the  distinction  between  open  and  closed  cracks  is 


determined  by  the  inequality 


which  can  be  evaluated  with  the  help  of  (28). 

The  results  presented  here  are  valid  only  if  this  inequality  is 
satisfied.  If  the  sign  is  reversed,  or  if  the  two  sides  are  equal,  then 
cracks  are  closed  and  do  not  affect  the  mechanical  and  thermal  response 
of  the  composite.  An  exception  arises  when  the  closed  crack  faces  slide 
in  shear.  Treatment  of  this  case  is  beyond  the  scope  of  the  present 
paper. 

Finally,  we  note  that  evaluation  of  stress  averages  In  the  phases 
is  trivial.  Since  open  cracks  do  not  support  any  stress. 

?0  1  7  *  2C  5  Q  •  (30) 
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2.4  Stiffness  at  Large  3 


When  a  fibrous  lamina  is  embedded  in  a  laminated  structure,  it  may 
crack  rather  extensively.  Under  such  ci rcumstances,  certain  components 
of  1.^  in  (3)  become  independent  of  6,  or  equal  to  zero.  It  is  worth 
evaluating  these  limiting  stiffnesses  since  they  correspond  to  the  worst 
possible  damage  to  a  particular  ply  of  a  laminate.  A  good  approximation 
to  the  actual  residual  stiffness  can  be  obtained  as 


l. .  a  lim  L. . 


(31) 


Now  as  s—  it  is  easy  to  see  that 


22’  44*  66 


whereas  the  remaining  components  of  the  compliance  matrix  remain  finite 
(and  equal  to  the  uncracked  value).  It  then  follows  from  simple  matrix 
inversion  that 


Hi  =  M33/y  *  L?1  "  ^2/L° 

★ 


h2  *« 


12'  '  22 


'22 


L23  =  0 


h*3  *  -M?3/y*L13  -  L?2LL/'-: 


(32) 


13 


'12  23'  22 


^*33  ■  <1^  *  -  (L,°,)2/L,° 


71 


33 


*23'  '  22 


-19- 


where  y  *  -  (^3)  •  Also, 


L55  *  L55 


L66  =  0 


For  a  moderately  cracked  ply,  say  3  ■  1,  a  reasonable  approximation 
to  the  stiffnesses  (ij  «  1,2,3)  is  still  given  by  (32),  trtiile  lg5  is 

of  course  equal  to  Lg5.  However,  and  Lgg  must  be  calculated  either 
directly  from  (3),  or  indirectly  from  (8),  (10),  and  (11). 


2.5  Selected  Results 


To  illustrate  the  evaluation  of  compliances  of  a  cracked  composite, 
and  of  the  "matrix"  strain  averages,  we  consider  a  T300  Gr-Ep  system. 
Table  I  lists  the  constituent  properties  of  fiber  and  matrix,  and  of  the 
uncracked  composite.  The  composite  compliances  were  obtained  from 
self-consistent  estimates  of  moduli. 

Figure  3  shows  changes  in  the  three  compliances  M22,  m^4,  and  Mgg, 
calculated  from  Equations  (8)  -  (11)  for  given  values  of  crack  density  s. 
Of  course,  all  these  components  increase  with  increasing  8,  but  their 
change  is  quite  gradual,  especially  at  low  6.  This  contrasts  but  is  not 
in  conflict  with  the  relatively  rapid  reduction  in  stiffnesses  which  we 
found  in  [5]. 

We  note  that  the  composite  without  cracks  is  transversely 
isotropic,  and  has  five  independent  compliances,  M.^.  When  cracks  are 
introduced  the  coefficient  M22,  M^,  and  Mgg  change,  the  material  is  no 
longer  transversely  isotropic,  and  the  number  of  independent  elastic 
compliances  increases  to  eight. 

The  changes  caused  by  s  in  components  of  strain  concentration 
factor  A  in  (22)  are  shown  in  Figure  4.  As  required  by  (26),  five  of  the 
eight  components  of  A  change  with  s.  Note  that  the  fiber  vo-lume  fraction 
appears  to  have  a  small  effect  on  A^.  However,  crack  density  s  can  have 
a  significant  effect  on  the  strain  concentration  factor  components, 
especially  at  relatively  low  values  of  S. 


3.  PROGRESSIVE  MATRIX  CRACKING  IN  A  PLY 


3.1  Preliminaries 


Consider  a  unidirectionally  reinforced  composite  ply  of  thickness 
2a,  which  is  embedded  in  a  laminated  plate,  shell,  or  a  similar  composite 
composite  structure.  The  structure  is  subjected  to  a  certain 
incremental  loading  program  which,  at  each  loading  step,  causes  a  known 
instantaneous  overall  strain  state  7  in  the  ply  under  consideration.  The 
strains  e^j  are  assumed  to  be  uniform  within  the  ply.  An  element  of  the 
composite  structure  which  contains  the  ply  is  shown  in  Fig.  5.  A  typical 
ply  is  bonded  to  laminated  layers  of  thickness  b'  and  b"  such  that  at 
least  one  of  the  adjacent  layers  is  much  thicker  than  the  ply,  if  b'  >  b" 
>  0,  then  b'  >  2a. 

At  a  certain  magnitude  of  applied  strain  F,  which  is  usually  much 
smaller  than  the  failure  strain  of  the  fiber,  a  system  of  matrix  cracks 
starts  to  develop  in  the  ply.  The  geometry  of  matrix  cracks  is 
influenced  by  the  strength  anisotropy  of  the  fibrous  ply,  and  by  the 
state  of  stress.  Clearly,  cracks  may  propagate  most  easily  on  planes 
which  are  parallel  to  the  fiber  direction.  An  unbroken  ply  segment  in  a 
typical  part  of  a  laminated  plate  or  shell,  which  is  removed  from  free 
edges  and  other  stress  concentrations,  usually  supports  a  state  of  plane 
stress.  In  the  coordinates  of  Figure  5,  this  suggests  that  components 
a22>  <7 23 »  and  °33  are  large  in  comparison  to  Matrix  cracks  parallel 

to  fibers  are  not  affected  by  033,  although  this  stress  component  may 
cause  cracking  on  planes  perpendicular  to  fibers.  The  remaining  stresses 


favor  cracking  on  x^-planes.  The  shear  components  a^2  ar>d  ®13  niay  be 
present  as  well,  for  example  in  bent  plates  where  their  magnitude  varies 
through  plate  thickness.  Again,  0^3  does  not  influence  axial  cracks 
parallel  to  the  fiber;  contributes  to  crack  growth  on  x^xj-planes, 
even  though  an  exceptionally  large  a ^  might  cause  cracking  on  planes 
inclined  to  x^  axis. 

These  considerations  suggest  that  the  matrix  crack  system  consists 
of  aligned  slit  cracks  of  width  2a,  which  are  infinitely  long  .in  the 
fi ber  di recti  on  X3. 

The  average  spacing  of  cracks  in  the  system  can  be  conveniently 
described  in  terms  of  the  crack  density  parameter  0,  the  average 
distance  between  cracks  is  2a/0,  Figure  5. 

Of  course,  the  actual  spacing  of  cracks  at  a  given  value  of  0  need 
not  be  regular.  Available  experimental  observations  suggest  that 
individual  cracks  are  initiated  at  randomly  distributed  locations  in  the 
ply,  and  then  propagate,  sometimes  in  an  unstable  manner,  across  the  ply 
width  [2,4].  The  number  of  these  cracks  grows  progressively  with 
increasing  7,  up  to  a  certain  saturation,  density,  typically  0  _<  1. 

Under  sustained  or  cyclic  loading  the  crack  density  may  increase  with 
time,  even  at  a  constant  stress  amplitude. 

Regardless  of  the  underlying  fracture  mechanism,  the  random  crack 
pattern  corresponding  to  a  given  average  density  0  can  be  illustrated  by 
the  sequence  of  Figure  6.  A  ply  segment  of  length  1  -  12.5  mm  and 


thickness  2a  =  0.125  mm  is  shown.  The  width  of  the  segment  is  assumed  to 
be  much  larger  than  2a.  In  agreement  with  experimental  observations,  we 
assume  that  the  actual  distance  between  any  two  adjacent  cracks  cannot  be 
smaller  than  the  ply  thickness  2a  [2-4].  Thus,  the  12.5  mm  long  segment 
can  contain  at  most  100  cracks  at  8  =  1.  To  illustrate  the  progressive 
cracking  process,  the  possible  crack  locations  in  the  ply  segment  at  8 
=  1  were  fixed  at  regular  intervals,  2a  apart.  Then,  a  sequence  of  100 
random  numbers  was  generated  and  each  number  was  assigned  to  one  of  the 
prospective  locations.  Four  values  of  8  were  selected,  and  at  each  8 
the  required  number  of  cracks  was  drawn  in  the  locations  with  the  lowest 
random  numbers.  For  example,  at  8  »  0.1  there  are  ten  cracks  in  the 
segment.  These  are  retained  in  their  positions  and  supplemented  by  20 
additional  cracks  at  8  =  0.3,  etc. 

The  sequence  shown  in  Figure  6  was  constructed,  under  the 
assumption  that  formation  of  a  new  crack  is  not  influenced  by  the 
presence  of  an  existing  crack.  Due  to  stress  relaxation  next  to  an 
existing  crack,  this  may  not  be  the  case  in  an  actual  ply.  At  low 
values  of  8  cracks  may  form  preferentially  at  locations  removed  from 
existing  cracks,  and  the  crack  distribution  in  the  early  stages  of  the 
process  may  become  somewhat  more  uniform  than  that  suggested  by  Figure 
6. 
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3.2  Crack  Interaction 


Before  the  onset  of  cracking,  the  local  strain  and  stress  fields  in 
a  typical  segment  of  a  thin  ply,  removed  from  concentrated  forces  and 
free  edges,  are  usually  uniform.  If  the  current  applied  strain  T,  and 
•uniform  thermal  change  0  are  known,  then  the  corresponding  ply  stress 
7  follows  from  constitutive  equations  (13). 

Once  a  particular  ply  begins  to  crack,  the  uniform  strain  and 
stress  fields  are  perturbed.  There  are  stress  singularities  at  crack 
tips,  and  the  stress  components  °13»  and  °23  vanish  on  each  crack 
face.  Nevertheless,  the  crack  affected  regions  are  localized  in  the 
neighborhood  of  each  crack.  Thus  in  the  early  stages  of  cracking,  each 
crack  can  be  regarded  as  essentially  isolated.  In  other  words,  the 
elastic  field  in  the  neighborhood  of  a  particular  crack  can  be 
determined  by  considering  that  crack  alone  in  the  ply. 

As  the  crack  density  increases,  so  does  the  totality  of  perturbed 
regions.  Ultimately,  the  cracks  will  interact  with  each  other,  and  the 
crack-affected  regions  will  eventually  occupy  the  entire  volume  of  the 
ply. 


It  is  important  to  estimate  the  distance  between  cracks  at  which 
crack  interaction  becomes  significant.  One  way  of  doing  this  is  to 
consider  a  ply  with  an  average  crack  density  S,  and  investigate  the 
differences  between  compliances  M  obtained  for  interacting  and 
noninteracting  cracks.  It  is  easy  to  anticipate  that  crack  interaction 
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will  cause  major  deviation  from  linearity  in  the  M22  component  which 
changes  most  rapidly  with  crack  density  6;  c.f..  Figure  3. 


Specifically,  with  reference  to  Equation  (4),  the  compliance  for 
interacting  cracks  is: 


0  1 

M22  *  M22  +  f  w 0  A22 


For  noninteracting  cracks,  M22  is  given  by  the  linear  formula 


m22  =  m22  +  J  *8  a22 


0  0 

where  M22,  A 22  are  values  of  M22»  A22  at  6  *  0. 


The  difference  between  the  above  compliances  can  be  expressed  in 

l  I 

terms  of  (M22  -  M22)/M22*  The  following  presents  this  ratio  for  the 
graphite-epoxy  ply  considered  in  Table  I;  the  magnitudes  are  not 
dependent  on  fiber  volume  fraction  Cf  for  0.2  <  Cf  <  0.6. 


0.25 

0.30 

0.40 

0.50 

0.105 

0.151 

0.271 

0.432 

These  results,  which  apply  to  cracks  already  in  existence,  suggests 
that  significant  crack  interaction  starts  at  e  =  0.25,  when  the  compli¬ 
ance  deviation  from  linearity  becomes  larger  than  10%. 

The  concept  of  crack  interaction  needs  to  be  understood  from  the 
stand-point  of  new  crack  formation  in  a  volume  of  a  ply  which  already  has 
an  average  crack  density  0.  If  a  crack  which  is  about  to  form  is  found 
in  a  ligament  between  two  existing  cracks,  it  can  be  regarded  as 
noninteracting  only  if  the  existing  cracks  are  at  least  8  or  more  ply 
thicknesses  apart,  i.e.  when  $  £  0.125  before  the  crack  forms. 

Of  course,  crack  initiation  at  random  locations  leads  to  consider¬ 
able  variations  in  local  spacing.  As  suggested  by  Figure  6,  individual 
cracks  may  be  much  closer  than  ply  thicknesses  at  M  0.125.  Therefore, 
one  can  conclude  that  crack  interaction  starts  very  early  in  the  cracking 
process.  Indeed,  assume  that  crack  interaction  is  present  for  all  values 


In  agreement  with  the  conclusions  reached  in  Section  3.2,  we  regard 
the  progressive  cracking  process  as  a  sequence  of  events  in  which 
individual  cracks  are  initiated  at  random  locations  in  the  ply,  and 
propagate,  possibly  in  an  unstable  manner  across  the  entire  loaded  width 
of  the  ply.  Formation  of  a  single  crack  represents  an  elementary  step 
in  the  process.  To  account  for  interaction  between  cracks,  we  assume 
that  the  stress  which  causes  the  crack  to  form  at  a  given  average  crack 
density  0  in  the  ply  is  equal  to  the  volume  average  (13).  Explicit  form 
of  relevant  components  is: 

—  _  o  __  _  o  . 

a12  *  2L66  6 12  3  2l66  +  ®  L66  A66)“X 

a 22  *  Li2(e"ll  -  aye)  +  \-ZZ^2Z  -  <*T«)  +  1-23(^33  *  aAe)  (36) 

*  ^**12^11  -  aT8)  +  L22  (*22  '  aT9^  +  L23^33  '  «A9^  1-22^22  ^ 

—  —  o  _  _  o  . 

<*23  3  2L-44e23  3  2L44e23(1  +  6  1-44*22)  ^ 

Fracture  criteria  for  crack  formation  during  progressive  cracking 
are  apparently  not  to  be  found  in  the  literature.  However,  it  is 
reasonable  to  expect  that  the  process  will  be  essentially  repetitive, 
governed  by  a  criterion  similar  to  that  for  first  ply  failure. 

In  general,  each  slit  crack  in  the  ply  can  open  in  three  modes, 
each  of  which  corresponds  to  one  of  the  three  stresses  (36).  To  fix 
ideas,  consider  that  the  slit  crack  extends  in  the  direction  parallel  to 
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the  fiber  axis  X3,  Figure  5.  Then  <112  in  (36)  causes  Mode  III  failure, 
033  Mode  II,  and  022  Mode  I.  Mode  III  can  be  excluded  in  most 
applications,  on  the  grounds  that  the  fracture  surface  in  a  fibrous  ply 
contains  numerous  longitudinal  asperities  created  by  the  crack  as  it 
finds  its  way  between  densely  packed  fibers.  These  asperities  are 
interlocked  and  thus  prevent  relative  shear  displacement  of  fracture 
surfaces  in  the  x^  direction.  Of  course,  Mode  III  and  ^12  are 
implicitly  absent  in  plane  stress. 

Modes  I  and  II  are  usually  interrelated,  but  to  a  different  extent 
in  different  composite  systems.  A  ply  failure  criterion  which  reflects 
this  fact  way  presented  by  Hahn  et  al .  [8]  for  the  case  of  >  °> 

°23  ^  0  in  the  form: 

(l-g)(T22/Tl)  +  9(T22/Tl)2  +  (723/yll)2  =  1  (37) 

where  g  is  an  empirical  interaction  parameter,  0  <  g  £  1.  For  g  =  1 
there  is  no  interaction  and  (37)  is  reduced  to  the  familiar  quadratic 
form.  The  and  are  critical  stresses  at  which  the  ply  actually 
fails.  These  stresses  increase  with  decreasing  ply  thickness.  Also,  as 
6  grows,  one  expects  that  formation  of  new  cracks  may  become  gradually 
more  difficult.  This  can  be  reflected  by  demanding  that  7j,  be 
nondecreasing  functions  of  S.  Thus  we  write: 

°I  *  <*l(a,B),  qj  »  oji  (a,0)  (38) 
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While  (37)  is  valid  for  o' 22  >  0,  the  composite  ply  may  also  crack 
if  ^22  <.  0  and  ^33  *  0.  Under  such  circumstances  (37)  should  be 
replaced  by 


(a23/7IlJ 


2 


1, 


(38) 


with  an  =  an  (a,  S,  <?22)* 


The. critical  stress  magnitudes  in  (37)  and  (38)  need  to  be 
established  experimental ly.  So  far,  only  first  ply  failure  magnitudes 
of  strengths  (38)  have  appeared  in  the  literature,  c.f.,  [3,4]. 
Information  about  7^,  especially  under  compression  (39),  is  generally 
not  available. 

If  (12),  (14)  or  an  equivalent  criterion  is  accepted  for  progres¬ 
sive  cracking,  then  the  average  crack  density  S  can  be  evaluated  during 
incremental  deformation  of  the  ply.  Suppose  that  the  composite  ply  is 
subjected  to  current  uniform  strain  7,  and  that  the  corresponding  crack 
density  is  80  >  0.  A  strain  increment  d7,  and  a  uniform  thermal  change 
de  are  applied,  and  a  new  magnitude  of  6  is  sought. 

Let 


e=e+de,  9 '  =  9  +  d9  (40) 

be  the  new  values  of  strain  and  temperature.  The  corresponding  average 
stress  in  the  ply  is,  in  analogy  with  (13): 


(41) 


o'  =L‘  e'  +  9'  l '  , 

****** 

where  l'  =  L'(8),  V  =  (B)  are  overall  instantaneous  material  properties 
of  the  cracked  ply  at  the  as  yet  unknown  value  of 

6  s  $0  +  d8 ;  dB  >  0.  (42) 

To  determine  if  d7  and  d6  in  (40)  will  actually  cause  additional 
cracking,  it  is  necessary  to  evaluate  first  stresses  (41)  with  L'  * 
l(B0),  i'  «  i (b0) ,  i.e.,  with  material  properties  at  the  original  crack 
density  80.  If  these  stresses  make  the  left-hand  side  of  (37)  or  (39) 
larger  than  unity,  then  an  increment  dB  is  required  to  reduce  L‘ ,  l'  to 
the  level  allowed  by  the  selected  failure  criterion.  When  (37)  or  (39) 
are  called  upon,  it  is  also  necessary  to  adjust  aj  and  according 
to  the  current  value  of  B. 

The  increment  dB  caused  by  (40)  can  be  evaluated  by  a  simple 
iteration  scheme  in  which  s  is  adjusted  so  that  the  stresses  (41) 
satisfy  (27)  or  (29)  while  7'  and  9'  are  applied. 

Once  the  new  magnitude  of  6  has  been  found,  it  is  Dossible  to 
evaluate  all  components  of  the  current  stiffness  and  compliance  matrices 
L  and  M. 

We  recall  that  both  L  and  l  are  decreasing  functions  of  B,  thus  the 
ply  softens  under  increasing  strain,  but  it  remains  elastic  during 


-31- 


f 


unloading  and  reloading  to  the  current  stress  level,  Figure  7.  This 
aspect  of  the  procedure  permits  evaluation  of  0  under  incremental 
deformation  even  when  the  critical  stresses  and  an  in  (37)  are  taken 
as  independent  of  0.  This  may  be  advantageous  in  practical  applications 
when  and  <7n  may  be  known  only  at  first  ply  failure,  i.e.,  at  0  0. 

If  the  first  ply  failure  values  are  taken,  the  predicted  magnitude  of  0 
will  probably  exceed  the  actual  one.  Therefore,  the  stiffness  of  the 
ply  will  be  underestimated  and  again,  the  model  predictions  will  be  on 
the  safe  side  in  most  structural  applications. 
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3.4  Results 


3.4.1  Stresses  in  a  Deforming  Ply 

To  illustrate  typical  results  provided  by  (37),  consider  a  graphite- 
epoxy  ply  with  properties  given  in  Table  I.  The  ply  is  deformed  by 
prescribed  strains  £22  and  £33  in  several  radial  directions,  while 
stresses  are  adjusted  so  that  »  0^2  ~  ^13  a  333  *  0.  For  greater 
clarity  of  presentation,  the  ply  strengths  are  taken  as  constant: 

-  _£r  _  _cr 

o\  a  °22  3  28.44  MPa,  ojj  *  023  *  37.95  MPa 

These  values  were  measured  for  first  ply  failure  of  an  unconstrained 
ply,  hence  they  represent  relatively  low  strengths. 

The  results  appear  in  Figure  8.  In  each  case,  the  applied  strain 
is  plotted  on  the  horizontal  axis,  the  average  stresses  and  crack 
densities  on  the  vertical.  The  constant  strengths  levels  are  plotted  as 
well.  It  is  easy  to  see  that  if  strengths  were  increasing  with  0,  then 
the  stress  components  would  be  elevated  accordingly. 

A  remarkable  feature  of  these  results  is  the  strong  effect  of  shear 
strain  £23  on  crack  density  increase.  When  only  shear  £33,  or  only 
transverse  normal  strain  £22  is  applied,  one  finds  that  6  increases 
about  twice  as  fast  in  shear  when  is  2e23  compared  in  magnitude  to  £32* 

At  the  same  time,  the  shear  strength  exceeds  the  transverse  tension 
strength  by  33%.  This  behavior  can  be  related  to  the  rate  of  change  of 
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the  relevant  stiffnessed  L22  and  L44  with  6.  We  recall  from  Figure  2  in 
[5]  that  L22  decreases  with  6  much  more  rapidly  than  L44,  from 
comparable  initial  levels.  Thus  023- is  relatively  large,  much  larger 
than  <122  for  2e23  >  £22*  ^ven  at  l°w  e23/£22  ratios  023  keeps  increasing 
with  applied  strain  while  always  decreases  after  onset  of  cracking. 

In  addition  to  stress  changes,  the  ply  experiences  stiffness  change 
under  deformation  conditions  of  Figure  8.  These  are  summarized  in  Figure 
9. 

3.4.2  Progressive  Cracking  in  Laminated  Plates 

The  above  procedure  for  analysis  of  crack  accumulation  in  a 
strained  ply  can  be  combined  with  laminated  plate  theory  and  thus  ex¬ 
tended  to  cracking  laminated  plates.  Figure  10  illustrates  stress  and 
stiffness  changes  in  a  plate  subjected  to  combined  loads. 

As  in  the  case  of  a  lamina,  cracking  in  laminated  plates  must  be 
analyzed  in  an  incremental  way.  In  many  applications  one  cannot  predict 
the  exact  loading  history,  or  a  detailed  analysis  may  be  undesirably 
complex.  Under  such  circumstances,  it  is  sufficient  to  estimate  the 
maximum  stiffness  loss  caused  by  the  greatest  possible  damage  in 
individual  plies.  Ply  properties  for  this  case  were  found,  for  large 
values  of  6,  Section  2.4. 


The  effect  of  a  fixed,  uniform  crack  density  in  all  plies  of  the 
laminated  plate,  on  longitudinal  (E|_)  and  transverse  (Ey)  in-plane 
moduli  of  plates  of  different  layups  is  illustrated  in  Tables  II,  and 
III  and  IV  [9].  It  is  seen  that  the  stiffness  loss  caused  by  extensive 
cracking  can  be  very  significant  in  certain  composite  systems. 

A  comparison  of  the  theoretical  results  with  experimental 
observations  is  presented  in  Figure  11.  Axial  stiffness  of  a  (0/903)5 
E-glass-epoxy  plate  is  plotted  against  observed  crack  density  for  a 
uniaxial  tension  test.  The  experimental  points  and  shear  lag  analysis 
were  presented  in  [10].  The  theoretical  (SCM)  curve  was  calculated  with 
the  self-consistent  approximation  described  earlier.  A  very  good 
agreement  of  this  curve  with  experiments  is  indicated.  A  remarkable 
feature  of  this  correlation  is  that  the  only  material  constants  used  in 
the  analysis  were  the  original  elastic  moduli  of  the  uncracked 
composite.  With  this  information  the  stiffness  changes  were  evaluated 
for  given  crack  density  from  (3)  for  the  90°  plies,  and  overall  laminate 
stiffness  change  was  found  from  laminated  plate  theory. 

Additional  comparisons  were  made  with  experimental  results  obtained 
by  Hwang  and  Hahn  [3]  who  measured  crack  density  as  a  function  of 
applied  stress  or  strain  in  AS/3501-5A  laminated  plate  coupons  in 
tension,  and,  with  similar  experiments  reported  by  Ryder  and  Crossman 
[11]  for  laminated  T300/5208  plate  coupons. 


Analytical  predictions  of  crack  densities  in  these  laminated  plate 
coupons  were  made  as  follows:  First,  laminated  plate  theory  was  used  to 
calculate  residual  thermal  stresses  in  each  ply  caused  by  cooldown  from 
the  curing  temperature.  The  temperature  difference  was  aT=147°C  for 
both  systems.  Then  ply  stresses  caused  by  a  uniaxial  applied  load  were 
evaluated  and  superimposed  on  the  thermal  residual  stresses  in  each  ply 
Next,  (37)  was  rewritten  for  each  ply  with  stresses  "a 22 ,  o*23  taken  in 
local  ply  coordinates,  X3  is  parallel  to  local  fiber  direction.  The 
°I»  aII  denote  critical  stresses  at  first  ply  failure.  The  parameter, 
which  accounts  for  failure  mode  interaction  in  an  approximate  way  was 
taken  as  g=0.1. 

Now,  to  determine  and  for  the  laminate  test  data  shown  in 
Figures  12  and  13,  one  calculates  the  local  stress  ^2  at  first  ply 
failure  in  the  90°  ply.  Since  033  =  0  in  this  ply,  one  obtains  <?i  from 
(37).  Then,  additional  stress  is  applied  *0  the  laminate,  and  crack 
density  3  in  the  90°  ply  is  calculated  incrementally,  according  to  the 
sequence  (37)  to  (42).  The  stiffness  of  the  90°  ply  gradually  decreases 
in  the  process.  When  the  applied  stress  reaches  the  magnitude  required 
for  first  ply  failure  in  the  +45°  ply,  ajj  in  this  ply  is  again  found 
from  (37). 

These  results  are  listed  in  Table  V  as  predictions  (A)  and  (B)  for 
the  two  coupons. 
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In  a  similar  way,  the  results  in  Figures  14  and  15  were  interpreted 
and  provided  values  of  and  7pi  n  predictions  (C),  (D)  and  (E) ,  Table 
V. 


The  experimentally  found  values  of  the  critical  stress  a j  and  op 
at  first  ply  failure  in  (37)  are  actually  the  sole  material  parameters 
needed  to  fit  the  data.  The  only  additional  information  entering  the 
analysis  are  original  elastic  moduli  of  the  undamaged  laminate,  and  the 
empirical  factor  g=0.1  in  (37). 

The  magitudes  of  and  are  fairly  consistent.  Note  that  pre¬ 
dictions  (A)  and  (8)  apply  to  both  0°  and  45°  plies  in  Figures  12  and 

13.  Differences  between  (A)  and  (B)  are  caused  by  differences  in  ex¬ 
perimental  data.  Prediction  (C)  is  sufficient  for  all  plies  in  Figure 

14,  except  for  the  -45°  ply  as  explained  below.  Predictions  (D)  and  (E) 
in  Figure  (15)  are  different  because  the  data  points  for  the  two  coupons 
are  so  far  apart. 

It  would  be  desirable  to  find  that  one  set  of  critical  stresses 
fits  all  data  for  a  particular  system,  or  for  the  two  grafi te-epoxy 
systems  of  Figures  12  to  15.  This  is  not  exactly  possible  but  we  note 
that  of  in  (A)  (B)  and  (D)  are  similar;  ideally  they  should  coincide 
since  they  refer  to  a  single  system.  Still  they  are  closer  together 
than  those  in  (D)  and  (E)  found  for  not  only  one  system,  but  also  one 
layup.  In  fact  the  two  coupons  A1  and  A2  were  taken  from  the  same 
laminated  plate. 
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Therefore,  one  can  conclude  that  the  differences  in  effective  ply 
strength  at  first  ply  failure  in  Table  V  are  caused  by  variations  in 
material  properties  that  affect  the  experimental  data. 

A  surprising  feature  of  the  experimental  results  shown  in  Figures 
12  and  14  is  that  no  cracks  were  observed  in  the  interior  -45°  plies. 
One  would  certainly  expect  to  find  cracks  in  these  layers  which 
experience  the  same  average  strain  state  as  the  +45°  plies  where 
numerous  cracks  were  found.  In  fact,  the  -45°  plies  form  a  single  ply 
of  double  thickness,  and  this  should  make  it  more  susceptible  to  matrix 
cracking  than  the  +45°  plies  of  single  thickness.  This  phenomenon 
requires  further  study. 

We  note  in  passing  that  not  all  experimental  data  can  be  approxi¬ 
mated  with  the  average  stress  condition  (37).  This  suggests  a  need  for 
a  more  accurate  method  of  evaluation  of  local  stresses  in  (37). 


3.5  Conclusions 


1.  An  analytical  technique  has  been  developed  for  modeling  of  progres¬ 
sive  transverse  matrix  cracking  in  laminated  composite  plates.  The 
analysis  consists  of  three  steps.  First,  self-consistent  estimates 
of  laminate  stiffness  changes  are  found  for  a  given  crack  density 
from  (3)  and  (4).  Next,  the  rate  of  crack  density  increase  under 
applied  load  is  evaluated  from  (37)  to  (42)  on  the  basis  of  a 
selected  ply  failure  criterion.  Finally,  the  ply  analysis  is  made 
in  terms  of  ply  stress  state  for  a  ply  embedded  in  a  laminate  with 
several  cracking  plies.  The  instantaneous  ply  stresses  follow  from 
a  simple  adaptation  of  laminated  plate  theory. 

2.  A  good  agreement  of  the  theory  with  several  sets  of  experimental 
data  was  found.  The  theoretical  predictions  require  a  minimum 
amount  of  experimentally  derived  information,  such  as  first  ply 
failure  stresses  in  the  cracking  ply,  and  elastic  moduli  of  the 
undamaged  laminate.  No  empirical  parameters  are  required  outside 
the  ply  failure  criterion  (37). 

3.  The  theory  makes  it  possible  to  calculate  the  instantaneous  crack 
density  in  each  ply,  the  instantaneous  stresses  in  each  ply  of  the 
laminate,  as  well  as  the  average  fiber  and  matrix  stresses  in  each 
ply,  after  application  of  each  load  or  strain  increment  to  the 
laminate.  The  laminate  can  be  subjected  to  combined  three 
dimensional  loading,  and  also  to  varying  uniform  thermal  changes  in 
the  course  of  mechanical  loading. 
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The  analysis  can  be  implemented  through  a  simple  numerical  routine, 
even  on  a  microcomputer.  In  fact,  the  entire  procedure  for  each 
loading  step  is  rather  similar  to  that  used  in  elasticity  analysis 
of  laminated  plates. 
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Table  I*  Constituent  Properties  and  Compliances  of  the  T300 
Gr-Ep  System 


Unit 

E33 

G31 

V31 

En 

G!2 

Symmetry 

Fiber 

(T300) 

1 03ksi 

33.00 

3.36 

0.410 

2.25 

0.78 

Transversely 

Isotropic 

103MPa 

227.5 

23.2 

15.5 

5.4 

Matrix 

(Epoxy) 

103ksi 

0.50 

0.19 

0.350 

0.50 

0.19 

Isotropic 

1 03MPa 

3.4 

1.3 

3.4 

1.3 

Composite  Compliances: 

Compliance  c^*0.2 

Cf*0.4 

cf=0.6 

I 


t' 

*?, 

0.2069 

0.1595 

0.1192 

P 

*1 

-0.1037 

-0.0779 

-0.0561 

M° 

13 

-0.0075 

-0.0041 

-  0.0028 

M22 

0.2069 

0.1595 

0.1192 

^23 

-0.0075 

-  0.0041 

-  0.0028 

M33 

0.0207 

0.0107 

0.0073 

r 

0.5108 

0.2726 

0.1244 

«S5 

0.5108 

0.2726 

0.1244 

mmm 

V  . 

M66 

0.6211 

0.4746 

0.3506 

r 

All  values  are  in 

units  of  (103 

MPa)"1 

L 
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*01  T 300/ 5208 
X  Of  WICRACXED  MOOULUS  EL 


less 

wm 

MM 

i  ,0 

(04,  90) 

100 

99 

99 

(02,  90) 

99 

98 

97 

(0  .  90) 

98 

96 

95 

(0  .  902) 

97 

92 

90 

(0  .  904) 

95 

86 

82  i 

(Oj»t*5) 

99 

97 

95  | 

(0  ,H5) 

98 

94 

92 

(0.(;4S)2) 

97 

91 

86 

(♦4S) 

91 

74 

59 

(02,;75) 

99 

96 

95 

(02,^0) 

99 

97 

95 

(02.*30) 

99 

97 

95 

(02.ri5) 

100 

99 

99 

401  T 300/5208 
X  Of  UNCRACXEO  MODULUS  ET 


lay-uiT''" 

MM 

m 

1.0 

(o2,:45) 

97 

91 

87 

(o,  :4S) 

96 

39 

84 

(0, ( *45 ) g ) 

95 

87 

30 

<:4S) 

91 

74 

59 

j  (02.;75) 

100 

99 

99 

(02,*60) 

99 

97 

96 

(02.*30) 

92 

79 

71 

(02,*15) 

38 

69 

59 

Table  II.  Theoretical  predictions  of  relative  changes  in  axial 
(El)  and  transverse  (ET)  Young's  moduli  of  laminated 
plates  at  given  values  of  crack  density  3  in  all  plies. 
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40*  boron/alumin ium 
X  OF  UNCRACKED  MODULUS  E, 


Lay-up''""''"' - 

0.2 

0.4 

0.6 

(04,90) 

96 

93 

91 

O 

o> 

#• 

CM 

O 

93 

37 

83 

(0.  90) 

89 

80 

73 

(0.  902) 

as 

71 

62 

(0.  904) 

80 

64 

51 

(04,:45) 

96 

93 

90 

(o2,:4S) 

94 

89 

85 

(0,  *45) 

91 

84 

78 

(0,(HS)2) 

89 

30 

73 

40t  BORON/ALUMIN ium 
X  OF  UNCRACXED  MODULUS  E, 


0.2 

0.4 

0.6 

(04,*45) 

95 

90 

87 

(02,;45) 

93 

36 

81  1 

(0.  *45) 

90 

32 

75 

(0,(:45)2) 

88 

78 

70 

Table  IV.  Theoretical  predictions  of  relative  changes  in  axial 
(El)  and  transverse  (Ex)  Young's  moduli  of  laminated 
plates  at  given  values  of  crack  density  3  in  all  plies 


PREDICTION 

al 

aIII 

MATERIAL 

LAYUP 

COUPON 

MPa 

MPa 

A 

64.5 

73.5 

AS/3501-5A 

(02/90/±45)s 

A-b-4 

B 

64.5 

125 

A-b-6 

C 

83 

105 

T300/5208 

(0/90/±45)s 

D 

45 

100 

(02/904)s 

A2 

E 

70 

100 

A1 

I 

^  Table  V.  Effective  Ply  Strength  At  First  Ply  Failure 


THREE-PHASE  MODEL  OF  A  CRACKED  LAMINA 


Figure  2. 

2A.  An  infinite  medium  with  aligned  fibers  and  slit  cracks, 
2B,  A  fiber  monolayer  with  cracks. 
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4. 

Changes 
"matri x" 

in  strain  concentration  factor  components  in 
caused  by  cracks  of  density  g  (T300  Gr-Ep). 

composite 

Components  (MPa) 


'.V' 


in  the  graphite-epoxy  ply  of  Figure  3,  as  functions  of  applied  strain. 


Stress  (MPs) 
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90*  plies 
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(0j/90/±45)e 
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833/10! 


*»"  0 
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O33  CMP*) 
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CMP*) 


033  (MP») 


Figure  10.  Calculated  changes  in  ply  stresses  5j j  and  plate 
stiffnesses  L-jj  in  a  graphite-epoxy  plate  loaded 
in  simple  tension  033  in  the  0°  direction.  All 
stresses  are  in  local  ply  coordinates,  the  fiber 
direction  coincides  with  local  X3  axis  of  a  ply. 


CO,  90  ,  90,90  3s  E  -  GLASS 


Figure  11.  Comparison  of  experimentally  observed  stiffness  changes 
and  crack  densities  [11]  with  theoretical  predictions. 
The  (SCM)  curve  was  calculated  from  the  self-consistent 
estimate  of  stiffness  in  90°  plies.  Equation  (3). 
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Figure  12.  Comparison  of  Predicted  and  Observed  Average  Crack 
Spacing  in  a  Laminated  Plate  Subjected  to  Uniaxial 
Tension.  Experimental  Data  from  Hwang  and  Hahn  (1982) 
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Figure  13. 


Comparison  of  Predicted  and  Observed  Average  Crack 
Spacing  in  a  Laminated  Plate  Subjected  to  Uniaxial 
Tension.  Experimental  Data  from  Hwang  and  Hahn  (1982). 
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